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We study the magnetic-field-induced quantum phase transition from a gapped quantum phase 
that has no magnetic long-range order into a gapless phase in the spin-1/2 ladder compound bis(2,3- 
dimethylpyridinium) tetrabromocuprate (DIMPY). At temperatures below about 1 K, the specific 
heat in the gapless phase attains an asymptotic linear temperature dependence, characteristic of 
a Tomonaga-Luttinger liquid. Inelastic neutron scattering and the specific heat measurements in 
both phases are in good agreement with theoretical calculations, demonstrating that DIMPY is the 
first model material for an S = 1/2 two-leg spin ladder in the strong-leg regime. 

PACS numbers: 75.10.Jm 75.40.Gb 75.50.Ee 



Gapped ground states comprising singlet pairs of spins 
are the prevalent nonmagnetic quantum disordered states 
in a variety of antiferromagnetic Heisenberg models 0- 
Among those models, two-leg spin-1/2 ladders with 
antiferromagnetic rung and leg exchanges, J rung and Ji eg , 
are the simplest whose ground states are yet non-trivial. 
These states give way to a Tomonaga-Luttinger liquid 
(TLL) — a critical state with fractional 5=1/2 spinon 
excitations — at a magnetic-ficld-drivcn quantum critical 
point (QCP) H. 

Although the quantum phase transition at such a 
QCP has been extensively investigated theoretically 0- 
8], there have been few experimental studies because 
of the scarcity of real systems with right energy scales. 
(CsH^N^CuCU, which was originally thought to be a 
ladder material [9[, later turned out to be a frustrated 
three-dimensional antiferromagnet [l(3]. In IPA-CUCI3 
111 Il2| . long range magnetic order — also known as a 
Bose-Einstein condensation of magnons [H, [H[ — due to 
interladdcr interactions dominates the magnetic-field re- 
gion above the QCP. Thus far, the only detailed re- 
port of a TLL in a two-leg spin-1/2 ladder has con- 
cerned (CRH 12 N) 2 CuBr4, a strong-rung material with 
Jieg/ Jrung ~ 0.25 [lj,[lal. For deeper understanding of 
ladders, development of new materials with a wide range 
of Ji e g/ Jrung will be crucial. Of special interest are mate- 
rials in the strong-leg regime, Jieg/ Jrung > 1; since quan- 
tum fluctuations are more prominent in this regime and 
as a result the singlets will be less localized, a state rem- 
iniscent of the resonating valence bond liquid 16|, Il7 |. 



In this Letter, we investigate a magnetic-field-induced 
quantum phase transition in (C7HioN)2CuBr4, DIMPY 
for short, a new material in which the CuBr^ 2 radi- 
cals form two-leg spin ladders along the crystallographic 
Our inelastic neutron scattering (INS) 
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a axis 

demonstrates that this compound is a spin-gapped quan- 
tum magnet with excellent one-dimensionality. Our 
specific-heat measurements reveal the presence of a 
TLL phase above the critical field H c = 3.0(3) T, with 
no long-range order at least down to 150 mK. With 
the aid of perturbative continuous unitary transfor- 
mations (PCUT) and state-of-the-art density-matrix 
renormalization-group (DMRG) calculations, we deter- 
mine the strengths of the rung and leg exchanges from 
the INS results in the gapped phase and the specific-heat 
results in the TLL phase with remarkable consistency, 
confirming that DIMPY is an ideal 5=1/2 spin-ladder 
system in the strong-leg regime. 

Single crystals of deuterated DIMPY were grown 
according to the method described in Ref. 
Prompt-gamma neutron activation analysis measure- 
ments showed that 67% of hydrogen sites are occupied 
by deuterium. The zero field INS experiment was per- 
formed on SPINS at NIST with a single crystal of a 3.5 g 
mass and a 0.5° mosaic spread. The measurements were 
made in the (h,k,0) and {hfl,l) reciprocal-lattice planes 
with a standard helium cryostat. The high-field INS ex- 
periment was performed on RITA II at SINQ, PSI. The 
sample consisted of two single crystals with a total mass 
of 2 g coaligned within 0.6°. The sample was oriented 
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FIG. 1: (color online). Dispersion measured by INS in 
DIMPY at T= 1.5 K as a function of ft, A;, and I. Lines in 
(a) and (b) indicate the gap energy. Lines in (c) are from 
PCUT calculations for an AFH two-leg spin ladder for differ- 
ent Values Of X = Jleg/ /rung- 



with the (h,0,l) plane horizontal and was cooled in a 13.5 
T vertical-field cryomagnet. The data rate was increased 
by employing a multi-blade crystal analyzer and a po- 
sition sensitive detector [lj|. A Be (or BeO) filter was 
placed after the sample to remove high-order contamina- 
tion, selecting a final neutron energy of 5.0 (or 3.7) meV. 
The specific heat measurements were made with relax- 
ation calorimctry at the NHMFL, Tallahasse, on a single 
crystal of an 8.2 mg mass in fields up to 18 T applied 
parallel to the c axis. 

Figure [T] summarizes the zero-field dispersion mea- 
sured at T=1.5 K by INS along three high symme- 
try directions in the reciprocal space [2(J. We per- 
formed a global fit of all collected data to a dynamic spin 
correlation function with the approximate spin-gap dis- 
persion t{q)=\J A 2 + v 2 sin 2 [27r(0.5 — h)] [2ll| . convolved 
with instrumental resolution, finding A = 0.32(2) meV, 
v = 2.36(4) meV. The individual data points shown in 
the figure were obtained by fitting a resolution-corrected 
lineshape to each constant- q (or constant-energy) scan. 
Note that Figs.QJa) andQJb) are shown on a much finer 
scale than Fig. [TJc). Within a scale as small as 20 /J.eV, 
dispersion is absent along the c direction, and only a 
very weak dispersion, if any, of at most 50 fieV is found 
along the b direction 22[, indicating that DIMPY is an 
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FIG. 2: (color online), (a) Background-subtracted constant- 
q scan in DIMPY at the magnetic zone center (0.5,0,0.9) at 
T=1.5 K for magnetic fields H=Q, 5, 8, and 13.5 T. (b) 
Background-subtracted constant hw — 0.7 meV scans along 
the (ft,0,1.7- 1.44ft) direction at T = 1.5 K and H = and 5 
T. The dotted line is a guide for the eye. In both frames, solid 
lines are fits to a dynamic spin correlation function with the 
approximate spin-gap dispersion relation [U, convolved with 
instrumental resolution, and dashed lines indicate zero. 



excellent one-dimensional (ID) system. 

We have calculated the dispersion of an 5=1/2 antifer- 
romagnetic Heisenberg (AFH) spin-ladder system, using 
PCUT [23| around the limit of isolated rungs. The scries 
in x = J\e S / ./rung is obtained in the thermodynamic limit 
24 1 and is extrapolated in terms of an internal parameter 
25 using Pade resummation, yielding reliable results for 
large x especially for q close to the magnetic zone center. 
The lines in Fig. [He) are the dispersion for different val- 
ues of x, calculated in conjunction with the accurate gap 
value A = 0.32(2) meV. Best agreement with the data is 
obtained for x = 2.2(2), indicating that DIMPY is in the 
strong-leg regime. 

Figure shows the background-subtracted 

constant-q scan at the magnetic zone center (0.5,0,0.9) 
at T = 1.5 K in different fields. The background 
was determined at zero field by making energy scans at 
q = (0.35,0,0.9) and (0.65,0,0.9), away from the magnetic 
zone center, with the same instrument configuration 
and by fitting the results to a Gaussian profile over the 



3 




0.0 




O v V o 

o vV Q o 
> |( ^ ooo 



o 


5T 


V 


8T 


o 


18T 



(b) 



0.0 0.5 1.0 1.5 2.0 2.5 
T(K) 

FIG. 3: (color online). Magnetic specific heat C m of DIMPY 
as a function of temperature T for (a) H < H c and (b) H > H c . 
In the latter region, the data have been plotted as Cm/T after 
subtracting the nuclear-quadrupole contribution (see endnote 
[27j )• vf is extracted from data between 0.3 K and the upper 
limit of the T-linear region, indicated by an arrow. Inset (1): 
Semilog plot of the H < H c data against 1/T. Lines are fits to 
Eq. [T] Inset (2): Field dependence of the spin gap obtained 
from the data. 



range where no magnetic excitation is present. At zero 
field, the resolution-limited peak indicates the location 
of the spin gap. Such a peak is absent at and above 5 T, 
indicating that the magnetic field drives the system into 
a gapless critical phase. 

To examine the magnetic excitation spectra at zero 
field and in the gapless phase, constant-energy scans were 
performed at T=1.5 K for huj = 0.7 meV as shown in 
Fig. [2jb), where a constant background term has been 
subtracted. These measurements were made along the 
(A,0, 1.7— 1.44/i) direction to maximize the structure fac- 
tor. The q resolution-limited peaks at zero field are from 
one-particle excitations. The low-energy feature in the 
gapless phase, at 5 T, is clearly much broader than the 
experimental resolution, suggesting that it arises from a 
two-spinon continuum, not from one-particle excitations. 

To augment the INS results, we measured the spe- 
cific heat at T < 2.5 K, as shown in Fig. [3] The phonon 
contribution was determined from the zero-field entropy 
S = J (C/T)dT and has been subtracted from the data 
at all fields. The nuclear-spin contribution has also been 



subtracted through a simultaneous fit to the data for all 
fields at temperatures below 700 mK. 

At zero field and 2 T, exponentially activated behavior 
is found, as shown in the first inset to Fig. [3ja) , providing 
additional clear evidence for a spin gap below a critical 
field. The specific heat of a gapped ID AFH quantum 
magnet in the low-temperature limit is given by j2(| : 
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where h is the number of gapped low-energy modes and R 
the gas constant. Fitted at fc^Tc A to this expression, 
the zero-field data yields A = 0.32(1) meV — excellent 
agreement with the INS result — and h/v = 1.26(2) [27l |. 
Taking v = 2.36(4) meV from INS, wc find n = 3.0(1), 
which unambiguously indicates the threefold degener- 
acy expected for a two-leg spin ladder. The field de- 
pendence of A is shown in the second inset to Fig. [3] 
a linear fit gives H c = 3.0(3) T in good agreement with 
A/(ffMs) = 2.8(2) T, assuming .9 = 2.0. 

Above H c , the specific heat shows remarkable behav- 
ior. There is no A-like peak, indicative of a phase tran- 
sition, at temperatures down to 150 mK and magnetic 
fields up to 18 T. FigureEJb) shows the magnetic specific 
heat divided by temperature, C m /T, at 5, 8, and 18 T. 
As temperature decreases, C m reaches an asymptotic T- 
linear limit, characteristic of TLL, before an upward de- 
viation sets in — probably a precursor of long-range order- 
ing due to weak intcrladdcr interactions [281 ]. The low- 
temperature specific heat of TLL is given by conformal 
field theory as 29L 30j|: 



C(T) = 7Ti?^77y 
3 Vf[H) 



(2) 



where vp, the Fermi velocity, is the velocity of the gapless 
excitations. Using this equation, we extract vf — 2.79(8), 
3.27(11), and 2.89(9) meV respectively from the specific 
heat at 5, 8, and 18 T. 

From these vf and A, we now determine x = 
Jieg/Jrung and J run g- First, we perform a density-matrix- 
renormalization-group (DMRG) calculations for 5 = 1/2 
AFH two- leg ladders [31|, in conjunction with finite-size 
scaling, and obtain vp/Jieg as a function of qjinH/ J\ rtr 

From 



and A/ J rung as a function of x 
and A = 0.32(2) meV from the zero-field 



for fixed x [32 
this A/J rung 

specific heat and INS, we find J\ eg — which is x Jrung — for 
each x. With these Jieg, we then normalize the exper- 
imental values of vf and plot them with the theoreti- 
cal results, as shown for x = 2 and 2.5 in Fig. 2] Fi- 
nally, comparison of experiment and theory in this plot 
yields x = 2.3(2), for which A/ J mng =0.409(6) and thus 
J rung = 0.78(6) meV. 

To summarize, DIMPY undergoes a quantum phase 
transition at H c = 3.0(3) T from a gapped phase to 
a Tomonaga-Luttingcr liquid (TLL). Inelastic neutron 
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FIG. 4: (color online). Field dependence of the Fermi velocity 
vf scaled with J\ eg . Filled symbols are DMRG results for 
x = 2.0 and 2.5. Open symbols are from specific heat, scaled 
with Ji g = 1.56 meV (x = 2.0) and Ji eg = 1.95 meV (x = 2.5), 
where Ji eg for each x has been found from calculated A/J rung 
(see the inset) and the spin gap A = 0.32 meV. For x from 2.0 
to 2.5, scaled experimental data lie in the shaded region. 



scattering reveals the excellent one-dimensionality of this 
material and provides a firm value of the spin gap, 
A = 0.32(2) meV, as does the specific heat. In the TLL 
phase, the specific heat attains characteristic T-linear be- 
havior, yielding the Fermi velocity of the gapless ex- 
citations for the first time in any laboratory TLL. We 
obtain J rung = 0.78(6) meV from A and vf, and the ex- 
change ratio x = 2.2(2) from the zero-field dispersion and 
2.3(2) from A and Vf- These are consistent with previ- 



ous estimates, J rl 



: 0.75 meV and x = 1.94, from mag- 



netic susceptibility [18| . Three independent experiments 
yielding the exchange constants with consistency and in 
excellent agreement with theory establish DIMPY unam- 
biguously as the first ideal realization of an 5=1/2 AFH 
two-leg ladder in the strong-leg regime, thus opening up 
an avenue for investigating the properties of such a ladder 
in this poorly explored regime. 
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